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Abstract. Homogeneous, nearly-isotropic Bianchi cosmological models are consid- 
ered. Their time evolution is expressed as a complete set of non-interacting linear 
modes on top of a Friedmann-Robertson- Walker background model. This connects the 
extensive literature on Bianchi models with the more commonly-adopted perturbation 
approach to general relativistic cosmological evolution. Expressions for the relevant 
metric perturbations in familiar coordinate systems can be extracted straightforwardly. 
Amongst other possibilities, this allows for future analysis of anisotropic matter sources 
in a more general geometry than usually attempted. 

We discuss the geometric mechanisms by which maximal symmetry is broken in the 
context of these models, shedding light on the origin of different Bianchi types. When 
all relevant length-scales are super-horizon, the simplest Bianchi I models emerge (in 
which anisotropic quantities appear parallel transported). 

Finally we highlight the existence of arbitrarily long near-isotropic epochs in models 
of general Bianchi type (including those without an exact isotropic limit). 
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1. Introduction 

Mathematically, the simplest models of the Universe are those which possess a high 
degree of symmetry. Spacelike slices of the Friedmann-Robertson- Walker (FRW) 
solutions possess the maximal six Killing vector fields (KVFs), reducing the number 
of degrees of freedom to just one, the curvature radius. A historical method of relaxing 
this assumption to study more complex models is to remove three KVFs corresponding 
to rotational invariance, leaving only a single group which acts simply transitively on the 
spacelike sections. This yields the set of "Bianchi" models, named after the classification 
scheme for three-parameter Lie groups PQ. One of the advantages of this approach is 
that the evolution (Einstein) equations reduce to ordinary differential equations. This 
means that the full, non-linear behaviour of Bianchi universes is open to study by 
analytic methods (e.g. dynamical systems approaches |2J) or numerical integration. 

Observational interest in these anisotropic models has receded somewhat since the 
publication of full-sky maps of the temperature anisotropies in the cosmic microwave 
background (CMB) by COBE [3] and WMAP [HE]. It is known [6] that minimal 
assumption^] coupled with the observed near-isotropy of the CMB suggest that our 
universe has been highly isotropic since the last scattering surface at a redshift of 
z ~ 1100. Of course, this does not preclude the existence of mild anisotropy at a 
level below (or bordering on) our current detection threshold. 

In fact known 'anomalies' in the CMB data have been shown to be mimicked by 
a weak signal corresponding to the pattern expected in VII^ universes [7] [S] [H]. Such 
models are, in the consensus view, unrealistic in that they require an unfeasibly large 
value of Qk i n light of strong constraints from other observations including the Gaussian 
random fluctuations in the CMB itself [ID] . Further, the expected polarisation signals 
from a Bianchi component are apparently incompatible with the limits on B-mode 
signals and TI?-correlations in WMAP [11] . On the other hand, the dynamical model 
employed in these studies ignored certain degrees of freedom in the interest of simplicity. 
One result is that the shear principal axes have a fixed alignment relative to the residual 
symmetry axis of the models - this restricts the appearance of the Bianchi component of 
the CMB to a form which is not entirely generic [TTJ H21 H3] • In a full analysis, one further 
finds that certain sub-models have oscillatory dynamics, rather than the monotonically 
decaying shear seen in restricted solutions. This certainly permits the polarisation 
constraint to be evaded [12] and may allow for a self-consistent anisotropic Bianchi 
universe to provide an improved fit to the CMB data over concordance models. However, 
a proper statistical analysis to support such an optimistic statement is currently lacking, 
and it is more plausible that, even with the additional freedoms, one will not find a model 
consistent with all known constraints. 

On the theoretical side, the Bianchi models remain useful pedagogical tools 
regardless of their observational status. Furthermore the emergence of anisotropic 

% It is worth noting that these assumptions include a temporal Copernican assumption; arbitrary 
large-scale anisotropies may be concealed from observation on a single time-slice. 
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pocket universes in string-inspired eternal inflation models [M] reminds us that the 
simple, globally isotropic picture of the present-day cosmos may yet be replaced. 

The work described here is intended to clarify the status of Bianchi models which 
are close to isotropy. Modern approaches to observational cosmology generally impose 
linear perturbations on a fully isotropic (FRW) background. Bianchi models which are 
almost isotropic must be expressible in this way: as small, 'homogeneous' (in a sense to 
be defined) perturbations on an appropriate FRW background. Such decompositions of 
Bianchi models into non-interacting modes has been made for individual situations in 
the past, and in some isolated cases can be made exact (e.g. Refs JT5l [16j [T7J [18j fl9]); 
however, here we exhibit a method for understanding all available homogeneous linear 
perturbations about an exactly isotropic background. The mode decomposition we 
present is very naturally suited to use in observational studies and has already been 
combined with the CMB transport equations in Ref. [H] to calculate consistently the 
full range of CMB patterns in nearly-isotropic Bianchi universes |12j . 

The rest of this work is structured as follows. In Section [2] we discuss the origin 
of the Bianchi model representations of maximally-symmetric spaces, highlighting how 
ambiguities arise in the notion of translations and the particular difficulty in resolving 
these ambiguities in an open space. We expand these ideas in concrete terms for 
each space in Section [3j with some further detail on alternative coordinate systems 
in |Appendix A In Section |4j we perform the decomposition of suitable Bianchi models 
into a background with linearised modes and discuss the classification of such modes 
alongside their relation to the automorphism formalism. We relate our results to 
previous work in Section [5j Finally, in Section [6j we outline the properties of a class of 
almost-isotropic Bianchi models which do not possess the fiducial isotropic limits. We 
provide a brief summary in Section [7j 

We use a metric signature — h ++ throughout. Latin indices represent spatial 
components (1 — > 3) while Greek indices represent spacetime components (0 — > 4). 
Indices relative to an ortho normal frame are given a hat (a, b, . . .). Units are adopted 
so that c = 1 = 811G throughout. 

The reader who wishes to skip the derivation and discussion will find a concise 
list of Bianchi perturbations in Table |2j their evolution and constraints described by 
equations (48) and (49), and coordinate expressions to expand such modes relative to a 
given coordinate system in Section [3j 



2. Homogeneous, anisotropic models 

2.1. Notational background; the traditional approach 

Given a manifold Ai with metric g, symmetries are mathematically expressed by the 
invariance of g under appropriate transformations of M.. Continuous transformations 
form a Lie group of operations on M; it is consequently possible to show that all 
continuous group operations are built from infinitesimal generators, and that invariance 
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under the infinitesimal generators is sufficient to prove invariance under the full group. 
Invariance of the metric g under infinitesimal transformations is thus encoded in the 
notion of a vanishing Lie derivative along the tangent vector field: 

C $ g = (1) 

where £ is a vector field everywhere tangent to the motion, expressing the movement 
of each point towards a neighbour. A vector field satisfying condition is known 
as a Killing vector field (KVF). The linearity of the Lie derivative means that linear 
combinations of KVFs are Killing. Furthermore, the Jacobi identities can be used to 
show that the commutator of two KVFs is also Killing, i.e. C^^g = 0. 

In the context of Bianchi models, homogeneity is defined by the existence of a closed 
algebra of three Killing vector fields which are everywhere linearly independent (so that 
the isometry group is said to act simply transitively on M). Intuitively, this expresses 
the existence of paths covering the entire space, along which one may walk from any 
point to any other while seeing no change in the metric. The closure of the algebra is 
expressed by demanding that 

&,fc] = -CS&, (2) 

where the Cfj are constants and we have introduced a minus sign for later convenience. 
For n = 3 spatial dimensions (which we will assume throughout), it is conventional to 
decompose the structure constants Cf, as 

C% = -cn5) + a$ + e tjl n lk , (3) 

where n 13 is symmetric. By reparametrizing the & with a linear transformation of 
the form & — > one may simultaneously require n = diag(ni, n 2 , n 3 ) and set 

a = (a, 0,0); in this scheme, the Jacobi identities reduce to an x = 0. The standard 
approach is to partition the space of allowed ni,n 2 ,ns and a into distinct regions or 
'Bianchi types'; models within each type are not necessarily isomorphic spaces but 
(up to parity) they are deformable into each other by continuous transformations of 
the metric at a point while keeping the Killing fields fixed. Conversely, models falling 
into different types are not related by such a continuous transformation; it follows that 
time evolution preserves the Bianchi type. We may further determine which of the types 
permit isotropic sub-cases, and these are listed in Table [TJ (See, for a concise derivation, 
Ref. plj.) 

2.2. An alternative approach: explicitly breaking maximal symmetry 

The explanations given above form the standard approach to defining the Bianchi spaces. 
However, for near-isotropic cases, it can be unclear how the symmetries thus derived 
relate to more familiar symmetries of the exactly-isotropic FRW cases: why does a given 
maximally-symmetric model arise within certain Bianchi classes and not others? 

To answer this question, it is instructive instead to start with the isometry group of 
the maximally-symmetric cases, and then inspect the origin of the Bianchi-homogeneity 
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Table 1. Bianchi groups with FRW limit and their structure constants in canonical 
form with right-handed parity for the KVFs. The final column is the comoving 3- 
curvature scalar in the (3 = FRW limit. 



and isotropy subgroups. This gives an alternative view of constructing near-isotropic 
Bianchi models: namely, we start with a maximally-symmetric 3-space, identify the 
three Killing fields defining homogeneity, and finally add a metric perturbation which 
is invariant under the action of these fields (but not under the remainder of the original 
isometry group). Dependent on the particular properties of these preferred Killing 
fields, the perturbed space then falls into a specific Bianchi class. (Alternatively, one 
can introduce perturbations which are invariant only under the isotropy subgroup, which 
would lead to a Lemaitre-Tolman-Bondi solution.) 

For a maximally-symmetric 3-space, a counting exercise shows there to be six KVFs 
of which precisely three must vanish at any single chosen point p. Since tracing the 
integral curves of these locally- vanishing KVFs leaves p fixed, the corresponding algebra 
elements must generate the isotropy group of that point. We can therefore immediately 
split the KVFs into two sets of three, {Rj} (rotations) and {Tj} (translations), where 
Rj| p = and Ti\ p ^ 0. Up to reparametrizations, the set of isotropy vector fields 
for p are uniquely defined. In fact, one can fix most of the reparametrization freedom 
by demanding (without loss of generality) that Rj satisfies the canonical commutation 
relations for rotation generators, 

[R;, Rj] = — ejjfcRfc. (4) 

On the other hand the translations {Tj} have much more freedom. While we can fix 
part of the reparametrization freedom by demanding (again without loss of generality) 
that the Tj form an orthonormal basis at p, there remains the freedom 

T, >£, T, • //• R,. (5) 

where p{ are constants. It is this freedom which ultimately leads to certain maximally- 
symmetric spaces falling into multiple Bianchi types (Table [TJ: there is no natural way 
to pick out a preferred set of generators £;|§] 

Not all choices of p\ are permitted: the three vector fields £i must be closed 
[expression ([2])] and also everywhere linearly independent. If the former condition is 

§ From a group-theoretic standpoint, the manifold is identified with the cosets of the rotation group; 
but because the rotations do not form a normal subgroup, there is no natural composition law for the 
cosets. 
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not obeyed, self-consistent (i.e. path-independent) transport will only be available for 
tensors which are at least partially rotationally symmetric, which is not our intention 
here (and leads either to full isotropy or the intermediate Kantowski-Sachs models). If 
the latter condition fails, the subgroup generated will not act transitively, which also 
implies rotational symmetry about some point. One may get further with determining 
the allowed commutation relations in a quite general framework, but in fact it is more 
enlightening to adopt explicit flat, open and closed models. These three types exhaust 
the possibilities; this can be seen, for instance, by considering the need for constant 
scalar curvature in maximally-symmetric spaces. We perform the explicit search for 
homogeneity groups in Section [3j finding exactly the groups of Table [T] from this 
alternative approach. 

Recalling that any tensor x will be defined as homogeneous if it is invariant under 
the action of the translational subgroup (C^x = 0), we can now ask how such a 
mathematical description compares with an intuitive notion of homogeneity. In flat 
space, for instance, physicists would almost automatically assume that homogeneity 
implies invariance under the Euclidean translational group (leading to a Type-I 
anisotropic model). The most natural way to extend this to curved spacetimes is to 
require our defining translational group to be generated by vector fields which are 
tangent to a congruence of geodesies. By enumeration, we will see that only Type-I 
(flat) and Type-IX (closed) models can have this property. 

Even where this can be arranged, anisotropic quantities on curved backgrounds 
(i.e. in the Type-IX case) cannot simply turn out to be parallel transported, since the 
latter arrangement is path-dependent. Nonetheless, there are two special properties 
of Type-IX and Type-I symmetries of maximally-symmetric models which make them 
attractive, as we now outline. 

First, for a Killing field to be geodesic, it is necessary and sufficient that it have 
constant norm i.e. d a (\$,i\ 2 ) = 0. Hence the translations described by such a field move 
each point by a constant geodesic distance and are known as Clifford translations. This 
property can be used to demonstrate the impossibility of constructing such Killing fields 
on an open background [20J. 

Second, the group of translations in the Type-I and IX cases are closed under 
rotation (the groups are normal within the full isometry group, with an ideal subalgebra: 
[Rj,£j] = —tijk$,k)- In fact, a set of transitive Killing fields on a maximally-symmetric 
space which are closed under rotation in this way necessarily describe a congruence 
of geodesies. The converse can also be shown to be true, i.e. a transitive set of 
Killing vector fields is necessarily closed under rotations if each field is tangent to 
geodesic congruencesjj[] An observational consequence is that the CMB in Type-I and 

|| On a maximally-symmetric 3-space with (constant) Ricci curvature R, one can show that V a £b = 
±^/R/6rj a bc£ c for a Killing field everywhere tangent to a congruence of geodesies. (Here rjabc is the 
alternating tensor.) Such a Killing field is therefore determined solely by its value at some point which 
we take to be p where the Ri vanish. By isotropy of the space, the Killing field that results from action 
of the Rj on £j is therefore equivalent to the field generated from the rotated version of £j at p, i.e. a 
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Type- IX models has only dipolar or quadrupolar temperature anisotropie^] [2U [TT] . 
(An exception to this statement arising in a locally isotropic limit falling outside the 
framework expounded here is explained in Section |6j) 

Conversely complex patterns in the CMB arise in Types V, VIIo and VII^ precisely 
from the lack of global rotational symmetry in the propagation equations for anisotropic 
quantities. We can now understand this to be a consequence of any of three distinct 
but equivalent statements: 

• the Killing fields adopted are not geodesic on the background; 

• the Killing fields adopted do not have constant norm on the background; and 

• the three-dimensional subgroup adopted to define homogeneity is not a normal 
subgroup of the full isometry group, i.e. it does not map onto itself under the 
action of rotations. 

To be clear, this lack of symmetry is in addition to the local anisotropies which will 
be introduced into the metric; it refers to the way such metric perturbations will be 
transported from a point through the space. In spacetimes formed from such models, 
observers armed with gyroscopes are able to detect directly the changing alignment of 
the Bianchi component of the CMB (or other anisotropic observables) as they move 
from location to location. Such effects can also be found in the VIIo flat model, and 
in some closed Type-IX models, but they are entirely unavoidable in the open Bianchi 
models. 



2. 3. Definition of the basis tetrad 

To formulate the dynamical description, we introduce the invariant triad e, obeying 

[e^-] = 0. (6) 

The e.i are not in general Killing but form a basis in which any homogeneous tensor, in 
particular the metric, will have constant components. The non- degeneracy of the basis 
at every point is guaranteed by the transitivity of the action. 

The ej are computed by selecting any basis for the tangent space at a point p and 
dragging out with equation We can choose to orthogonalize the at p (and hence 
everywhere) relative to the background metric. We can further satisfy ej| p = £i\ p by 
a suitable reparametrization of the £j. One then has [e^e^] = -t-C^e^. By making 
further linear reparametrizations of the and the are brought into canonical 
form without disturbing the orthogonality (see Table [T]). We have intentionally left the 
e$ un-normalized, so that, in the maximally-symmetric background, g(ei,ej) = e 2a Sij 
for some a. Alternatively, one could set a = to gain an orthonormal triad, but the 
structure constants cannot generally be chosen to be canonically normalized in this case. 

linear combination of the This establishes closure under rotations. 

% For the first order CMB the geodesies can be taken from the background FRW spacetime [21] . This 
being conformally related to K (g) M. , the photon null geodesies follow the pattern of spatial geodesies 
on M. 
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Four-dimensional spacetimes with maximally-symmetric 3-spaces can be con- 
structed in the standard manner, by introducing a foliation of the 3-spaces labelled 
by a time t, and using the hypersurface-orthogonal one-form h = dt to complete the 
spacetime metric: 

g = -n(g)n + e 2a (t) e i g, e 3§.. ; (7) 

where {e*} form the basis dual to {ej}. Both {e 4 } and {e^} are time-invariant, i.e. their 
Lie derivative along n - the vector corresponding ton - vanishes. The tetrad {n, e;} 
is known as the time-invariant tetrad for the Bianchi spacetime [22J. 



3. Killing fields and invariant vectors for specific spaces 

In this section, we write explicit expressions for all six Killing vector fields of each 
maximally-symmetric manifold, then perform the decomposition into isotropy and 
homogeneity groups as described in Section [2] 

This both elucidates the origin of the Bianchi symmetries and allows construction of 
coordinate expressions for the Killing fields and invariant basis. While such expressions 
have been calculated by Taub [23J (see also Ryan and Shepley [21]), they are presented 
in a form which results in a zero-order metric not immediately familiar to cosmologists. 
We will present our results here in as generally applicable a form as possible; longer 
expressions which result from writing the vector fields in familiar coordinate charts are 



given in |Appendix A 



3.1. Flat space: Bianchi Types I and VIIq 

The six Killing vector fields of flat space described by Cartesian coordinates {x 1 } are 
most simply expressed as the translations Tj = di and the rotations -R, = e^kX^d^, with 
commutation relations 



[r*,r y ] = o 

[Ri, Rj] = —tijkRk 

[Ri,Tj] = —eijkTk. (8) 

Solution of the closure requirement Q in terms of the p{ ^ leaves substantial freedom, 
but up to reparametrizations the two possibilities obtained give the Type-I and VIIo 
KVFs (£( and ^ Io respectively) in canonical form: 

e? = Tt (9) 

and 

£Vii =Tl + R 1 

= Tj j ± 1. (10) 
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Figure 1. The Killing fields of Type VIIo compared to those of Type I; some integral 
curves of both are shown (as thick and thin solid lines respectively). All the Type-I 
Killing fields (along with £^ n ° and £3 °) induce translations along the coordinate 
axes, whereas £^ n ° (as a linear combination of d\ and R\) induces a more complex 
spiralling motion. The vectors of the invariant triad e^ n °, e^ n ° are also plotted. These 
spiral in response to the motion of the Killing fields; see equation (111. 



The Tj have constant norm, form a subalgebra which is also closed under rotations and 
have geodesic integral curves; as noted in Section [2j these conditions are all equivalent. 
The Type- VIIo KVFs do not satisfy such conditions, because integral curves of the £1 
fields describe a spiralling motion (see Figure [lj). 

Because the background is flat, parallel transport of homogeneous quantities 
between different points on the manifold is uniquely defined; it may be verified that 
in the Type-I (but not VIIo) case the Lie transport coincides with parallel transport. 

In the Type-I Abelian and so form their own reciprocal group. For 

Type VII , the reciprocal group is given by solving [ej,£j] = subject to the constraint 

£-i|iE=0 £i|a;=0- 
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These solutions are illustrated by plotting the vectors e 2 and e^ n ° at various points in 
Figure [TJ Recalling that anisotropic quantities in the Type- VIIo case will have constant 
components relative to the e^ n ° basis, the origin of spiralling patterns in such universes 
becomes clear from inspecting this diagram. 
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Figure 2. The open space as a hyperboloid embedded in Minkowski space; here a 
cut is taken along a 3 = 0, and a portion of the resulting geometry illustrated in two 
projections. In particular we plot integral curves of KVFs, decomposed about the 
point p (where a = 1) as described in the text. The rotational motion is shown as 
a dotted line, while boosts are shown as thin solid lines. However, the boosts do not 
form a closed group and thus do not qualify as generators of homogeneity according 
to the criteria in Section [2] Thus the Bianchi models take a linear combination of the 
boosts and rotations to form the basis set of Bianchi KVFs at p (see equation 16 1. 
Here we have plotted integral curves of the Bianchi Type-V motions (heavy solid 
lines) which agree with the boosts in the tangent space of p, but deviate on scales of 
order the curvature radius. The rotational asymmetry of these integral curves, which 
is unavoidable in open models, ultimately leads to the non-trivial appearance of the 
Bianchi CMB. 



3.2. Open space: Bianchi Types V and Vllh 

To obtain coordinate expressions for the Killing and basis vectors, we consider the well- 
known embedding of an open manifold in Minkowski space. The ambient metric is 

ds 2 = -d(a ) 2 + dia 1 ) 2 + d(a 2 ) 2 + d(a 3 ) 2 . (12) 

The hypersurface representing the open space, with curvature set to 1 without loss of 
generality, is defined by 

1 = (a ) 2 - [(a 1 ) 2 + (a 2 ) 2 + (a 3 ) 2 ] . (13) 



The group of continuous global linear transformations which preserve (12 13) is the 
restricted Lorentz group SO + (3,l). There are six Killing vector fields which can be 
decomposed about a point p into three rotations and three boosts (Bi): 

Bi = a% + a°d h (14) 
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where the special point is taken without loss of generality to be a = 1, a 1 = (i > 0). 
These Killing fields have well-known commutators: 

[Ri,Rj] = — €ij k R k 
[Bi, Bj] = tijkRk 

[Ri, Bj] = — €ijkB k . (15) 



The commutation relations (15) show that the boosts Bi are not suitable generators 
for homogeneity because they do not form a closed group. One must therefore 
consider linear combinations of the form (|5]). It may be verified that the only closed 
possibilities (unique up to transformations under the isotropy subgroup {Ri}) give either 
a realisation of the canonical Type-V or VII^ KVFsj^] In the former case we have 

^ = B,, ^ = B 2 -R 3 , ^ = B 3 + R 2 . (16) 

The effect of this choice is shown in Figure [2] by displaying the hyperboloid described 
in the space spanned by the a 1 coordinates (taking a 3 = 0). In particular, we have 
illustrated integral curves through p of both the B%, B 2 subset (thin solid lines) and 
the £i, £ 2 subset (heavy solid lines). Both of these sets span the a 1 , a 2 tangent space at 
p which allows for a direct comparison. While the Bi subset is closed under rotations, 
this is not the case for the £j subset; for instance, it is clear from the figure that £2 is 
not obtained from any rotation about p of £1. We emphasize again that, despite their 
attractive qualities, the Bi are not suitable for transporting anisotropic quantities since 
they do not form a closed subalgebra (Section [2]). 
For Type Vllh one has 

Zj 1Ih = Vh(Y (j + 1). (17) 



Note the similarity of the relation (10) between Type-I and VIIo Killing fields with the 
relation (17) between Type V and VII^. Both introduce a 'spiralling' motion into one of 
the Killing fields. In the flat case, this is fixed (in the canonical decomposition) to have 
a unit comoving scale length which is later adjusted in physical scale by the choice of 
spatial scale-factor at redshift zero. In the open case, there is an additional freedom in 
the relative scale of the isotropic curvature and spiral motion; this is fixed by adjusting 
the value of the structure constant \fh. 

We can find invariant vectors by solving ej] = subject to the constraint ej = £j 
at a = 1. For Type V, we find 

= (e x - a ) d + (e x - a 1 ) 3 1 - a 2 3 2 - a 3 d 3 

= e % v , (18) 

+ For both Types V and Vllh, the normalization of £j for j 7^ 1 can be changed independently of the 
normalization of £1 while still leaving the structure constants in canonical form. In the VII^ case the 
scaling of £2 must be equal to that of £ 3 . Here we have chosen normalizations which fit in with our 
formalism, in that they lead to a spatial metric oc Sij in the exactly isotropic limit, cf. equation |7|. 
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where x = — ln(a° — a 1 ). (For a geometrical interpretation of x, see Appendix A.l 
Vllfc, a suitable invariant basis is 
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This simple formulation arises because the & are closed under the action of R±, so 
the Jacobi identities applied to = [Ri, [$,j,e k ]] may be used to show that the {ef} 
are closed under the action of R±. The corresponding commutation coefficients are 



fixed by the behaviour of Ri on the tangent space at the identity, 



i.e. 
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e^. (These may be verified directly from equation 18 ) The 



construction is therefore completely analogous to that of the VIIq invariant fields. 



3.3. Closed space: Bianchi type IX 

The closed case has been discussed extensively in Ref . [19] . A suitable embedding may 
be considered with the ambient metric and surface function 



ds 2 = (da ) 2 + (da 1 ) 2 + (da 2 ) 2 + (da 3 ) 2 
1 = (a ) 2 + (a 1 ) 2 + (a 2 ) 2 + (a 3 ) 2 , 



with KVFs 

Ri = ^ijkOpdk 

Ti = a l d - a°di. 
The commutation relations are 



(20) 



(21) 



[Ri, Rj] 
[Ti,Tj]-- 



-ijk-' 



\Ri, Ta 



■Rk 
'^ijk-Rk 

1 /J " ~^ijk^k- (22) 

In the vicinity of a = 1, a* = 0, the Tj correspond to Euclidean translations along the 
ith coordinate direction and the -Rj correspond to rotations about aP = 1, a 1 = 0. 
There are two simply-transitive subgroups 

$f = \{Ti± Ri) , (23) 
and both yield Type-IX structure constants: 

(24) 



With the factor 1/2 in (23 ), the £ 4 have canonical commutation relations, while reversing 
the sign of the makes this set canonical too. One may verify that both sets of 
Killing vector fields have constant norm on the surface, so are tangent to congruences 
of geodesies as described for the flat case above. From this it follows that parallel 
transport and Lie transport along a Killing field can differ only by a rotation in the 
plane perpendicular to the transport direction. This behaviour can succinctly account 
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for the appearance of the CMB in Type-IX models: the temperature pattern is a 
pure quadrupole but a screen-projected polarization vector spirals at a constant rate 
in conformal time |12j . (However, for an exception to this behaviour see the near-VIIo 
models in Section |6j) 

Since = 0, the invariant fields are also Killing (the metric is bi- invariant; it 

is, actually, the Killing metric for this space). For consistency we should adopt = £+, 
e- x = £r, although these roles are reversed under parity 



4. Enumerating the linearised Bianchi modes 

In this section, we present a derivation of the independent modes which could arise in 
an almost-isotropic Universe consistent with Bianchi-type symmetries. 

Our starting point is the metric in equation ([7]). We can break isotropy while 
preserving homogeneity by adding a spatial metric perturbation of the form 

Sg = 5g ab (t)e a ® e b . (25) 

As described in Section [2j the perturbed metric is invariant under the action of the 
homogeneity group, which is a subgroup of the original isometries for the unperturbed 
FRW spacetime. The full, perturbed metric can be written in the form [2"2l |2"5] 

g = -h®h + e 2a{t \e 2m ) ab e a ® e b , (26) 

where the symmetric, trace-free (STF) f3 a b{t) and, generally, a(t) are perturbed from 
their FRW values. By construction, the isotropic limit corresponds to (3 ab = (3 ab = 0, 
where an overdot denotes the derivative with respect to t. The evolution of /3 contains 
all information about the anisotropic perturbation (at any order). 

The evolution (Einstein) equations are most simply expressed relative to an 
orthonormal frame, denoted e„. One may construct such a frame explicitly |22j via 
the transformation 

ea = e-(e^) aa e tt . (27) 

The expansion H, shear cr &b and Fermi-relative rotation of the orthonormal frame^] 
are defined by 

e a ■ [n, e g ] = -H5 db - a hh - e hh £l\ (28) 

where <r a g is symmetric and trace-free. We shall now linearise about the background 
FRW solution by requiring 

= 0(e) 

P ab /H\ = 0(e) 

e- 2a \P ab C 2 /H 2 \ =0(e), (29) 

* Note that f.g = e.ggfi for comparison with certain works such as those by Hawking and 
Collins P3 123 [|>1. 
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where e< 1, The final condition can be interpreted as putting a lower bound on the 
wavelength of the perturbations as a fraction of the horizon size, which is necessary for 
spatial derivatives of perturbed quantities to remain small. We perform a fiducial fluid 
decomposition of the stress tensor as 

= pn^nv + p{g fW + n^nv) + 2n^P v) + ir^, (30) 

where p is the fluid energy density measured by observers comoving with the n 
congruence, p the isotropic pressure, P M the momentum density and the anisotropic 
stress. Here P^ and tc^ u are orthogonal to n (P^n^ = n^n 1 ' = 0) and tc^ u is symmetric 
and trace-free. We further assume 

\n, u /H 2 \ =Q(e), \PjH 2 \=0{e), (31) 

and will in practice shortly restrict = 0, i.e. a perfect fluid, although unlike 
many Bianchi model analyses we will allow for small non-zero momentum density P M 
(corresponding to tilt of the fluid source). 

The first-order kinematic and dynamical quantities in the orthonormal frame are 

e a aa =a & - (3a a a a + 0(e 2 ) (32) 

e a fi &h = n hi + (3^ + (3 ib n* b + 0(e 2 ) (33) 

°JH = P-JH + 0(e 2 ) (34) 

H =a (35) 

Q £ /H = 0{e 2 ). (36) 

These equations are non-tensorial (valid only in the orthonormal frame); the components 
of a and n on the right-hand sides are the canonical, time-invariant structure constants 
and the components of (5 ab are as defined in the time-invariant frame. 

The expansions above may be inserted directly into the orthonormal-frame Einstein 
equations (see, e.g., Ref. [2]; or for a somewhat different form, Ref. |26j), giving 

p = 3H 2 + 3 R/2 - A + 0(H 2 e 2 ) (37) 

H = -H 2 - (p + 3p)/6 + A/3 + 0(H 2 e 2 ) (38) 

Kh = -3#4s - % b + + 0(H 2 e 2 ) (39) 

Pa = e- a (3kba b - e a6c /W) + 0(e~ a HCe 2 ), (40) 
with the fluid conservation equation 

p = -3H{p + p) + 2e~ a a~ a P h + 0{e~ 2a HC 2 e 2 , H 3 e 2 ), (41) 

where the structure constants a and n again take their canonical numerical values. In 
equations (37) and (38), 3 R and 3 S db = 3 R &b — 3 R3 ab /3 are the isotropic and anisotropic 



parts of the spatial 3-curvature respectively: 

3 Re 2a = 3 R a fi e 2a = -n ab n ab + ^n aa n bb - 6a a a a 



+ (3 ac (-4n ab n cb + 2n ac n bb + 12a a a c ) + 0(C 2 e 2 ) (42) 
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Any STF (5 has S = 0, TZ 
Any STF (5 has S = 2, TZ 



0. All modes transverse. 
0. All modes transverse. 



Table 2. The decoupled linearized modes of Bianchi perturbations about an FRW 
Universe, expressed as the perturbation to the spatial metric relative to the {e a } basis. 
Under each mode, the eigenvalue for the anisotropic curvature (S) is given along with 
the perturbation to the isotropic curvature (TZ) and momentum-density constraint (P) 



as described by equations (46 1, (51) and (49) respectively. For transverse modes, the 
momentum density is zero. 



ab 



2a 



^ n c{a n b)c n cc^(aS) 



+ 2 



% e cd(a n b)c a d 



- e cd{ah)e U ^d ~ t cd{ a n b)e a dPec + ^ cd( a n b)c a ^ed + 2/5 ' cd n c(a ^ d 

+ 2 Pd(a n b)c n cd ~ Pcdn cd n {(li) - n cc n d{h (3~ b)d \ + 0(C 2 e 2 ), (43) 

where angle brackets (()) denote symmetrising and removing the trace over the enclosed 
indices. 

By construction the zero-order contribution to the curvature arises from the FRW 
background, with Ricci scalar values -Rfrw listed in Table [T] The anisotropic curvature 
to first order is linear in (5 so that one may write 

(44) 



3 S ±t e 2a - S cd& b(3cd, 3 Re 2a - 3 R FRW e 2a + 



ab 



and equation (39) becomes 

lb + afffltf + e- 2a s rMh p cd - K lh 



(45) 
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Anisotropic stresses could arise, for example, from collisionless free-streaming, magnetic 
fields or topological defects; see Ref. (27]. Their effect will be dependent on the source 
evolution equations for tt^; by combining such expressions with our present formalism, 
many existing analyses of anisotropic stresses could be extended from the simple Type-I 
setups usually employed in such studies. However, for the present work we will assume 
a perfect fluid (7r-g = 0); then the Bianchi perturbations decouple into independent 
eigenmodes satisfying 

°abcdPab — ° Pab °ac°bd \^°) 

where the are listed for each type in Table [2] The split into three classes (s, v 
and t) reflects a rotational symmetry of the structure constants about the e.\ axis. The 
modes transform irreducibly as scalars (s; spin-0), vectors (i>; spin-1) and tensors (t; 
spin-2) under this symmetry]!] these evolve independently in our linear approximation. 

The appearance of complex eigenmodes for Type-VII^ t modes indicates a time- 
dependent rotation of the mode expansion relative to the n-induced point identification 
between time slices. Real physical quantities are obtained by combining the complex 
conjugate mode pairs in linear combinations in the initial conditions; one may verify 
that the reality is conserved by the time evolution. 

A general Bianchi linear perturbation is written 

Pab = £ AW^W, a, h = AM(t)fi$\ (47) 

m=l m=l 

with the amplitude of each mode evolving independently as 

A (m) // + mA {m) 1 + S (m) A (m) = q ^ 

Here we have made use of conformal time r]; primes (') denote d/dr] = e a d/dt and 
H = He a = a'. Those modes with 5^0 describe damped harmonic oscillators and 
in general undergo oscillations, continually exchanging anisotropy between shear and 
curvature as they die away. The shear of other modes (S = 0) decays as e~ 3a ; these 
modes also admit an = constant solution which is a pure gauge artefact (see 



Section 4.1 ). 



The tilt constraint, equation (40), reads 



Pa = e- a J2A^P^ m \ (49) 



m=l 



where 



Pt ] = WSPat ~ WfrW - (50) 

The components are listed for each mode in Table 2 Those modes with vanishing 
momentum density are labelled transverse in the table. The projected divergences of 

jj This decomposition can be made in the non-linear case |28j and is an example of the automorphism 
approach to Bianchi dynamics (e.g. [29] and references therein). 



Linearization of homogeneous, nearly- isotropic cosmological models 



17 



the shear (i.e. the right-hand side of the tilt constraint) and the metric perturbation 
vanish in such modes at linear order. 



The Friedmann and Raychaudhuri equations (37) and (38) are modified by the 
perturbation to 3 R at first order, 

5 

5( 3 Re 2a ) = ^ {m) A (m) 



m=l 



=n ab pif, (si) 



with 



TZ ab = -4n ac n cb + 2n ab n cc + 12a a a b , (52) 
although in practice most of the IZ^ vanish (Table ^1). The resulting backreaction on 



the evolution equation (45) is second order and is therefore not considered. 



4-1. Gauge modes 



The perturbed metric in equation (26) is not fully gauge-fixed by the requirements we 
have outlined (that \(3 a b\ be small; that the commutators of the e a are canonical; and that 
[to, e a ] = 0). There are two residual gauge freedoms: a global time shift in the labelling 
of the homogeneous hypersurfaces, t — >■ t + 5t, and a constant linear transformation 
of the e a that preserves the canonical form of the structure constants. Applying these 
gauge freedoms to the background FRW metric generates apparent perturbations that 
are actually gauge artefacts. 

Small time shifts do not affect (3 ab at first order, but global reparametrizations 
e a — > T ab e b generally do. (The latter transformations correspond to time-independent 
automorphisms in the group-theoretic terminology.) Under reparametrization, the 
structure constants transform tensorially; hence a a and n ab transform as components 
of a covector and pseudo-tensor respectively: 

a a -> T ab a b (53) 
n ab -> det (T)(T~ 1 ) ca (T~ 1 ) db n cd . (54) 

Here, we are concerned with reparametrizations that maintain the canonical form of 
a a and n ab . For such a T ab close to the identity, we can write T ab = 5 ab + t ab where 
t ab = 0(e). This transformation generates a non-zero, constant (3 ab = t( ab ) to first order 
in e. To see how these gauge modes relate to the eigenmodes in Table |2j note that any 



mode with S = admits a solution of equation (46) giving constant {3 (hence vanishing 
shear) and, by construction, the anisotropic curvature is vanishing. Such a model is 
therefore isotropic, hence FRW, and the "perturbation" (3 ab is pure gauge. Working 
through the Bianchi types with FRW limits, we can enumerate the allowed t ab that 
preserve the canonical structure constants and hence determine the gauge modes. For 
example, in Type V the allowed t ab have t a ± = but such transformations can generate 
an arbitrary constant (STF) f3 ab . This is consistent with the findings in Table [2] that 



all (3 ab have S = in Type-V models. For the other four types, we similarly find that 
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the gauge modes correspond to linear combinations of the eigenmodes in Table [2] with 
S = 0. There are no gauge modes in Type IX. 

The scale-factor e a also generally changes under reparametrization of e a ; to linear 
order, e a — > e a [l + t aa /3}. In non-flat models, the change in scale-factor under a gauge 
transformation is compensated by a non-zero lZ a bPab to preserve the intrinsic curvature 
3 R of the homogeneous hypersurfaces. For the gauge modes in Types V and VII/j, 
t aa = — 3/3n and, since 3 R is unchanged, we require lZ a b(3 a b = 12/3n (Type V) or 
Ttabfiab = 12/i/3n (Type V). These are consistent with the values 1Z = 24 and 24h 
reported in Table [2] for those eigenmodes with gauge-mode counterparts and 1Z non- 
zero. 

4-2. Classification of modes 

We can classify modes according to their transformation under suitable rotations of the 
basis vector fields. Under e a — > R a b^b, with R a bR ac = 5b c , fiab transforms tensorially, 
i.e. (3 ab — > R a cPcdRbdi thus preserving (3 ab = in isotropic models. Any mode is either 
invariant under such a transformation, or transforms into an equivalent mode. 

For Types I and IX, the structure constants are invariant under any rotation. This 
means that all modes can be made from a single basis mode under a combination of 
rotations and linear superpositions; hence all modes must have the same dynamical 
behaviour. However in Types VII , V and VII^ the structure constants possess only a 
residual U(l) symmetry corresponding to rotation about e.%. In Table [2j the s modes are 
invariant under U(l), whereas v\ and are a pair of modes which transform into each 
other under a spin-1 representation. In the case of ti and t 2 for models other than VII^, 
the modes transform into each other under a spin-2 representation. For the Vllh spaces, 
the diagonalization procedure applied to the anisotropic curvature introduces complex 
elements, describing left- and right-circularly-polarized gravitational wave modes which 
do not mix but individually transform according to a complex spin-2 representation. 

Our classification of the Type-V and VII modes into s, v and t classes is analogous 
to the Fourier-based classification of linear perturbations about flat FRW models into 
scalar, vector and tensor perturbations, whereby scalar, vector and tensor modes 
transform as spin-0, 1 and 2 representations under U(l) rotations about the Fourier 
wavevector. More generally, cosmological perturbations over maximally-symmetric 3- 
spaces can be decomposed into scalar, vector and tensor modes following [30]. Adopting 
the synchronous gauge, the perturbed line element is 

ds 2 = -dt 2 + a 2 (t)[(l + 20) 7ij + h ij ]dx i dx j (55) 

for a background FRW spacetime with spatial metric a 2 (t)^ij in some coordinate chart 
{t,x z }. The trace-free decomposes as 

hij = (ViVj - 7 iiV fc V fc /3)£ + 2V (l B j} + (56) 

where Vj is the spatial covariant derivative compatible with 7^-, S is a scalar field, 
Bi is a solenoidal 3-vector field and Wy is a STF and transverse 3-tensor field. The 
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Mode Vector generator B 



Vvi 


e 2 






e 3 




VII vi 


( v^e 2 - e 3 ) / (h -\ 


-1) 


v 2 


(e 2 + Vhe 3 J 1 (h H 


-1) 



Table 3. The generators for the vector modes in Types V, VII and VII^ relative to 
the invariant basis. For the Type-VIIo generators, h should be set to zero. Note 
that in all types displayed, e.i and are solenoidal. The e\ direction does not feature 
because it is either Killing (Type VIIo), in which case Vrj-Bj) vanishes, or not solenoidal 
(Types V and VLI h ). 



perturbations (j> and 8 describe scalar perturbations, Bi vector perturbations and 
tensor perturbations. In the closed case, the decomposition into scalar, vector and tensor 
parts is unique but in the non-compact case this is only true if the perturbations satisfy 
appropriate asymptotic decay conditions. 

Our t-mode metric perturbations can all be shown to be transverse and so can be 
interpreted as tensor perturbations. However this interpretation is generally not unique 
since the required asymptotic decay conditions are not met here. Similarly, our f-mode 
metric perturbations can be interpreted as vector modes, i.e. they can be written as 
V(j-Bj) for a solenoidal vector B^. For the w-modes, it is always possible to take Bi to 
be homogeneous; explicit constructions are given in Table |3j 



The metric perturbation from (3^ (see equation 57) can be shown to be transverse 
for all Type- IX modes and, since the spatial sections are compact, these are uniquely 
to be interpreted as tensor perturbations (gravitational waves pj5]). One can also 
show (Section [3]) that the invariant fields are Killing, which explains the lack of any 
homogeneous vector perturbations. 

The ambiguity in the scalar, vector and tensor decomposition in the non-compact 
case is simply illustrated in Type I since the background is flat space. All modes are 
transverse but may nevertheless be expressed in terms of vector modes: consider, in a 
coordinate basis {x, y, z}, the Type-I mode with only flu = 02i — 1 non-zero. One may 
decompose = d^Vj) with V = (xd y + yd x ). The vector V is solenoidal (V • V = 0), 
but also curl free (V x V = 0), so that it too can be further decomposed into the 
scalar V = V0, <p = X V- Similar statements can be made for any Type-I mode. These 
decompositions are possible in spite of standard results because does not decay 
towards infinity. 

4-3. Metric perturbation in coordinate bases 

The expressions in Table [2] can be expanded relative to a coordinate basis by using 
the explicit expressions for the triad of invariant vector fields given in Section [3] and 
Appendix A 
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The Type-I case is trivial, since the invariant vectors are simply the coordinate 
differentials, ei — di. We therefore consider, for illustration, the next simplest case, that 
of the VIIo modes on a flat FRW background. In this case, expressions (11 ) can be used 
to generate an anisotropic synchronous-gauge metric perturbation, 

8 9lxv = 2e 2a p ab (e a ),(e b ) u , (57) 

for each decoupled linear Bianchi mode. 

For example, the VIIo V\ perturbation has an anisotropic correction to the 
metric with coordinate-basis components {bg\i,bg\z) proportional to (cosx 1 , sina; 1 ). 
For the VIIo s perturbation, the anisotropic metric perturbation has components 
(5<7u, 8g 2 2, Sg3s) proportional to (2, —1, —1). In this case, we can regard the mode also 
as a Type-I perturbation, i. e. the metric is invariant under the action of both the Type- 
I and VIIo groups and the perturbed space is locally rotationally symmetric, having 
four KVFs ({Tj} and Ri). Finally, the VIIo t\ perturbation produces an anisotropic 
metric perturbation with 5g 2 2 = —5g33 oc cos 2a; 1 and 5g 2 ^ oc sin 2a; 1 . Physically, this is a 
standing wave resulting from the superposition of two counter-propagating left-circularly 
polarized gravitational waves with wavenumber 2. 



5. Relation to non-linear solutions 



Compared with previous approximate approaches to analysing the behaviour of near- 
isotropic Bianchi models, the advantage in the description above is its use of the physical 
anisotropic curvature and shear variables as natural linearisation coordinates. The 
linearisation applies when anisotropic shear, stress and curvature are all small; this 
is expected to be suitable for describing any residual anisotropy in the Universe today 
but may break down at arbitrarily early or late times. Therefore we now place our 
evolution equations in the wider context of present knowledge of the properties of the 
exact solutions. 



5.1. Early-time attractors 

For a given shear at the present epoch, it is usually supposed that anisotropies were 
larger at sufficiently early times when the FRW scale-factor was smaller; ultimately as 
ol — y — oo one expects to enter a non-linear regime in which at least one of the conditions 



(29) breaks down. While this expectation is broadly confirmed, we will discuss below 
that there are a few exceptions. 

The supposition of a non-linear phase arises from the early-time attractor for open 
and flat models [3T] , which is a Kasner shear-dominated approach to the singularity. In 
that case, all isotropic and anisotropic curvature lengthscales become super-horizon and 
the Universe is indistinguishable from a Type-I model with shear diverging as (1 + z) 3 
for z — > oo. Our description does not explicitly resolve the Kasner vacuum behaviour 



because we have neglected a 2 terms in the Friedmann constraint (37) and Raychaudhuri 
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equation (38). The overall dynamics of such a model may be understood by matching 
an early Kasner phase onto our late-time linear description. 

However it need not be the case that lengthscales of interest are super-horizon at 
early times: for instance in Misner's closed 'mixmaster' (Type-IX) Universes [32] , quasi- 
Kasner phases are interspersed with 'bounces' which reverse the sign of the shear. These 
bounces occur when the universe becomes sufficiently small along two of its principal 
axes that the curvature re-enters the horizon. (Pancaking universes, in which only one 
axis becomes short, collapse to a singularity without a bounce; this corresponds to 
travelling along the + direction in Misner's terminology.) 

In fact even in open and flat models, the Kasner attractor may be avoided: Lukash 
[TT] showed that in a VIIo or VII^ model the initial conditions may be chosen so that the 
anisotropy in gravitational waves is hidden in super-horizon curvature at early times. 
In our notation this corresponds to setting A^' = 0, A® = Aq at an initial singularity. 
Inspecting equation (48) and its eigenvalues (Table [2]) shows that, with such initial 
conditions, shear is generated in the tensor modes only once the mode enters the horizon 
(i.e. r)y/S > 1), in agreement with Lukash's description. The dynamical behaviour 
is then correctly modelled all the way back to the initial singularity using our linear 
approximation. 

The existence of these modes has been highlighted by various authors, for instance 
Collins and Hawking (who refer to the existence of a 'growing mode' due to the behaviour 
at very early times [26]) and Barrow and Sonoda (Ref. |18J, Section 4.9; the regular t 
modes and decaying s modes are seen, but not the v modes because the tilt is set to 
zero). In our own work we have previously [12] denoted these modes t% (distinct from 
the shear-diverging to). Note that, contrary to the notation in that work, the present 
paper's t\ and t<i modes are rotations (or conjugates) of each other, both able to give 
regular high-redshift behaviour if the initial data specifies A 1 = 0. 

Similar tuning of the initial conditions can be applied to produce closed Type-IX 
models with a non-oscillatory approach to the initial singularity, despite their generic 
mixmaster behaviour described above. As with Type VII, the existence of regular Type- 
IX early-time solutions can be understood through the decomposition into an exact 
background with super-horizon gravitational waves; this separation has been further 
explored in Refs. [IE] and [T9~] . 



5.2. Late-time attractors 

We have considered above the 'early-time' behaviour, by which we mean the behaviour 
when isotropic and anisotropic curvature scales are super-horizon. We will now turn to 
the 'late-time' behaviour of models after the characteristic scales are sub-horizon. We 
should note, however, that the real Universe could be in either of these states - i.e. 
there is no guarantee that we have now entered a 'late time' in terms of the dynamical 
evolution of anisotropics in the Universe (if such anisotropies indeed exist). 

Near-isotropic Type-I, V and IX models are relatively uninteresting at late times. 
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Type-I and V models which are sufficiently close to isotropy will always tend to full 
isotropy at sufficiently late times [26J. This is reflected in the trivial behaviour of the 
linear dynamical modes, all of which have 5 = (Table [2]) and so have shear that 
decays as (1 + z) 3 . Since Type-IX models recollapse at late times, they do not achieve 
a late-time limit but instead re-enter a time-reversed 'early-time' phase. 

Of more interest are the Type- VII models since, according to Collins and Hawking 
[26J, the subset which ultimately isotropize is of measure zero in the initial conditions. 
This conclusion is reinforced by the more recent numerical analysis of Wainwright [33J ; 
the late-time behaviour of the Type- VII models therefore warrants closer inspection. 

The quantities which define isotropization for each mode are the expansion- 
normalized shear (£( m ) = A^ m ''/'H) and curvature (S^A^fH 2 ). Isotropization occurs 
if both of these tend to zercfpfl as t — > 00; conversely the universe becomes arbitrarily 
anisotropic (and our linear approximation breaks down) if either quantity diverges. We 
will first consider the behaviour of flat VII models, followed by the open VII/j models. 
In both cases only the t modes are of interest, since all other modes have no interaction 
with the anisotropic curvature and the shear decays as (1 + z) 3 . 



For t modes of type VII at late times, equation (48) shows (for 7 = p/p + 1 > 2/3) 
that oc e 2M? e a( - 37//2 ~ 2 ). In the case of radiation domination, 7 = 4/3, the expansion- 
normalized shear thus undergoes constant-amplitude oscillations; at second order the 
amplitude has been shown to decay logarithmically [31 J. Such slow decay has been 
emphasized by Barrow [3U [27] in the context of nucleosynthesis constraints, but it is 



important to realize that it applies only to the late-time expansion in which \J S^r) 1. 
This will be attained during radiation domination only if the current physical length 

1 /2 

scale of the spiralling is smaller than th/Qmi where r# is the present-day Hubble 
length and Qr and Qm are the current radiation and matter density parameters. If 
this condition is not satisfied, the shear and anisotropic curvature of the tensor modes 
were essentially decoupled during radiation domination because the scale-length was 
super-horizon. 

In Type-VII/i models the conformal Hubble parameter % is constant but non- vanish- 
ing at late times once isotropic curvature dominates. (The following considerations 
therefore do not apply to models with dark energy.) In our linear approximation the 
solutions for the t modes in the curvature-domination limit are 



A{rj) = exp <^ -U 



1 ± ^i-s^/h 2 

with W = +Vh and <S( m ) = 4(1 — iy/h). The amplitude of the oscillations is 
determined by 1 ± Re a/1 — S^fH, which is equal to 1 ± 1 for all h. Consequently the 
amplitude of the expansion-normalized shear oscillation is either constant or decays as 
e~ 2a . The 'frozen-in' constant-amplitude mode was previously noted in the analysis of 
Doroshkevich et al. 



ft As noted by various authors, this restriction is stronger than required by the observational constraints 
on isotropy; see the last paragraph of Section [6j 
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5.3. Implications of nucleosynthesis in the observed universe 

The above considerations have implications for constructing a Bianchi model which is 
compatible with big-bang nucleosynthesis (BBN) constraints. Any remaining difficulties 
in fully reconciling observed abundances with theory based on the fiducial FRW picture 
(e.g. Ref. [35]) are relatively minor; BBN thus constitutes significant evidence against 
any model which produces an expansion history significantly different from FRW at 
T ~ 10 9 K. In fully non-linear relativistic evolution, shear enhances deceleration, 



modifying the Friedmann (37) and Raychaudhuri (38) equations so that 



p = 2»H 2 - a 2 + 3 R/2 - A (59) 

H= -H 2 -\{p + ?> P ) + ±-\a\ (60) 
o 3 3 

where a 2 = o~ h yi ah /2. For this reason, significant shear at high redshifts is prohibited by 
our knowledge of the chemical composition of the Universe. Since for most linearized 
modes a oc (l + z) 3 (this relation also holding exactly for the non-linear Type-I models), 
even weak shear constraints at high redshift translate into stringent upper bounds at 
the present day. Accordingly, for many models, BBN is known to constitute a strong 
limit (competitive with or even stronger than microwave background bounds - see [36]). 
However, exceptions arise in certain models which, in terms of our description, either 
have regularized initial conditions (the t modes in which shear anisotropy vanishes and 
anisotropic curvature scales are super-horizon at high redshift), or much slower than 
(1 + z) 3 decay at late times (the logarithmically decaying modes for very small spiral 
scales). In such circumstances, the BBN constraints offer far weaker limits on present 
day shear; it is these models, therefore, which should be of most interest to CMB 
phenomenologists [12]. 



6. Other Bianchi types 

We conclude with some comments on the Bianchi models which do not possess a strictly 
isotropic limit (Types II, IV, VIo, Vlh and VIII), but nonetheless include models which 
are close to isotropy. By definition these models can still be thought of (inside the horizon 
for a limited epoch) as perturbations on a maximally-symmetric FRW background, but 
the symmetries of the perturbations will no longer exactly correspond to Killing fields 
of the background maximally-symmetric models. 

For near-isotropy on any single timeslice, the expansion-normalized shear cr^/H 
and anisotropic curvature 3 S^/H 2 must be 0(e). Suppose we now take a general 
Bianchi model, not necessarily with a strict FRW limit, which on some timeslice satisfies 
these requirements. On this timeslice we may, without producing a physical change, 
reparametrize the triad of invariant vectors so that (3 a b is zero. This may bring the 
commutation constants C c ab into a non-canonical form, but we can at least rotate the 
basis vectors so that n a t> is diagonal and dj = for % > 1. We then form a tetrad 
with the timelike normal and extend to nearby timeslices using the same time-invariant 
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prescription described in Section |4| At times sufficiently near our original chosen 



timeslice, the conditions described by (29) will apply so that our perturbative expansion 
of the dynamics is valid. 

The physical quantity 3 S & ~ b /H 2 expressed on an orthonormal tetrad remains 
invariant (up to rotations) under the above transformations. Since on the preferred 



timeslice we have set j3 a b — 0, only the zero-order part of expression (42) for this 
quantity contributes. It follows that the structure constants of the time-invariant tetrad 
must be equal to those of a 'nearby' isotropic Bianchi type (for which the zero-order 
contribution to 3 S~f ) would vanish), plus a small perturbation. In such a case of near 
isotropy, all those Bianchi types that do not have a strict isotropic limit must fall into 
one of the following three possibilities. 

(i) All of the structure constants are small in the sense that 

e- 2a \C 2 /H 2 \ = 0(e). (61) 

The nearby isotropic case is Bianchi Type I, and the actual universe could be any 
Bianchi type. Because the structure constants are small, the anisotropic part of the 
CMB [11] appears as a pure quadrupole at first order provided the near isotropy 
holds between recombination and the present epoch. The dynamics are modified 
such that 

«S + 2U ^al + ^ = °> ( 62 ) 
where <S s g = 3 S~f ) e 2a is time independent at first order. 

In those models with an exact FRW limit, we are describing a situation where 
the shear is small and the isotropic curvature scale and any twist scale are small 
compared to the Hubble radius. In other models, while the linearization applies, 
we can say the Universe is in a flat quasi-Friedmann epoch, in agreement with the 
terminology of Ref. [2] (see Section 15.2.1). S &b can be made arbitrarily small by 
appropriately scaling the structure constants. The epoch of apparent isotropy can 
thus be made arbitrarily long in any Bianchi type, although this implies fine-tuning 
of the initial conditions. 

(ii) a± = 1, ni — 0; n 2 and n% are both small in the sense that 

e- 2a ai \ ni /H 2 \ < 0(e) (63) 

The nearby isotropic case is Type V. One possibility for the actual Bianchi model 
under analysis is Type VII^ (signn2 = sign 713). Building a rescaled triad shows 
that this captures VII^ models in the limit where the twist scale is large compared 
to the Hubble radius. 

The other possible types which can achieve near-type-V isotropy are Type VI^ 
(signn2 = —sign 723) or IV (one of n2 or n% vanishes). In all types, the expression 
for the anisotropic curvature reads 

3 S, b e 2a = Sl^ cd + (n 2 - n 3 )/3* = (n 2 - n 3 )/^, (64) 
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which drives the corresponding gravitational wave term in Table [2] according to 

A (t»)H + 2UA {t2) ' = -(n 2 - n 3 ) . (65) 

The universe can then be said to be in an open quasi- Friedmann epoch, 
(hi) |n 2 — 1| and \ns — 1| are small; a\ = 0; ni is small such that 

e" 2 >2 - nl\/H 2 = 0(e), e~ 2a \n x {n 2 + n 3 )\/H 2 = 0(e). (66) 

The nearby isotropic case is Type VIIo and the actual Bianchi model under analysis 
is Type VIII {n x < 0) or IX (m > 0). 

Building a triad with canonical Type-IX commutation relations from a near- VIIo 
model will produce a decomposition with highly anisotropic metric (\(3 a b\ large), 
violating the original linearisation conditions (29). The near- VIIo models therefore 
form a distinct isotropic limit which only the alternative analysis of this section will 
correctly capture. 

For all near- VIIo types, the modified curvature expression reads 

* s «s - O*- " % (67) 

which drives the scalar mode, i.e. 

A {s) " + 2HA^ s)l = -m . (68) 
3 

(Note that if only this driven scalar mode is present, the values of n 2 and n.3 are 
unobservable since (3 a {, has a U(l) symmetry about the e x direction; the physical 
situation is then identical to a near-Type-I model with small n\.) 

We have already commented that in the near-Type-I case, the anisotropic-sourced 
part of the CMB remains a pure quadrupole if the linearisation holds all the way back 
to recombination. Given the constraints on flatness in today's Universe, we also know 
that the near-Type-V case could be seen only in a near-flat limit which, likewise, would 
involve only a quadrupole contribution to the CMB. On the other hand the near- VIIo 
models (case iii) give rise to more interesting insights, as described below. 



6. 1 . Near VIIo ' spiral patterns in Type-IX models 

Case (iii) of the preceding text shows that spiral patterns in the CMB can arise in 
Type-IX models. Specifically, by perturbing n\ about a VIIo background with t or v 
modes on some initial hypersurface, one will generate a closed Type-IX universe which 
has observable spiral structure in the CMB. Note that, because it is the difference 
\ n 2 ~ "^il 1 ^ 2 which must be small compared to the inverse comoving horizon scale He a , 
the VIIo spiral scale (set by the magnitude of e a /n 2 or e a /n 3 ) can be sub-horizon. 
This is a surprising conclusion: Grishchuk showed [16] that all Type-IX models can 
be decomposed into an exact FRW background and maximal-wavelength gravitational 
waves - the origin, in this picture, of an apparently independent spiral scale is at first 
unclear. 
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However, one can confirm the spiralling even in the Grishchuk decomposition by, 



for instance, transforming a model obeying restrictions (66) into a frame in which the 
structure constants are canonical. In this "Grishchuk frame", (3 a b is nonlinear even 
though the shear and anisotropic curvature remain small while the mode is on super- 
Hubble scales. The exact geodesic equation (e.g. equation 23 in Ref. [11J) can be 
applied; because the Grishchuk-frame components of e 2 ^ are far from the identity matrix, 
the geodesic behaviour is necessarily altered from the canonical isotropic Type- IX model. 
Once the resulting geodesic is expressed relative to an orthonormal group-invariant 
triad, one sees multiple spirals per Hubble length, in exact agreement with the effects 
derived from the original near-VIIo frame. In other words, the spiral scale information 
is correctly encoded in the non-linear gravitational waves of the Grishchuk picture. 

Before leaving this topic, we should note that spiralling in Type-IX has recently 
been discovered independently by a quite different approach [37] . 

6.2. Isotropic curvature and anisotropy 

For the near-V and near-VIIo models (cases ii and hi), the horizon re-entry of isotropic 
curvature will unavoidably trigger the appearance of anisotropy in the universe, because 
the anisotropic- and isotropic-curvature scales are geometrically linked. For instance, 
consider the Type-IX models which are passing through a near-VIIo phase. To first 



order, the Friedmann constraint (37) reads 

m 2 = pa 2 - ni . (69) 

Consequently a detection of the universe's closure is possible once the expansion has 
slowed sufficiently for ni/% 2 to exceed some small threshold. A similar criterion applies 



to the detectability of the driven scalar shear mode (68), if its magnitude is initially 
close to zero. So at the same time as the isotropic curvature enters the horizon, the 
scalar-mode driving term will become important. 

In the linear approximation, there are no late-time solutions with small anisotropics 
in any of the models introduced in this section, unless 7 < 2/3. At sufficiently late times, 
the near-isotropy will necessarily break down for these classes of models, unless some 
form of dark energy is dominant or the expansion is reversed. In terms of the exact 
Bianchi solutions for 7 > 2/3, the near-isotropic solutions form saddle points since they 
have both stable directions (decaying shear modes) and unstable directions (growing 
modes, those which are forced by the small 3-curvature). 

It has previously been emphasized [31], however, that the existence of long- 
lived quasi-FRW solutions in all Bianchi types is all that is required to make those 
Bianchi types of observational relevance; even without dark energy there would be no 
observational need for the isotropization to be asymptotic as t — > 00. All this makes it 
impossible to rule out the relevance of any given Bianchi type to the observed universe, 
although the observational near-flatness places partial constraints on the scale of the 
Bianchi geometry relative to today's horizon. 
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7. Summary 

We have investigated in some detail the geometrical origin and time-dependent 
behaviour of Bianchi perturbations on maximally-symmetric (i.e. FRW) backgrounds. 
Many results obtained relate closely to previously-known properties, but to our 
knowledge this work is the first to consider the various models in a unified way using 
physically transparent variables throughout. 

By diagonalising the linearised Einstein equations, we have enumerated the linear 
modes in such models, and have shown each to be a scalar, vector or tensor mode on a 
suitable background. It is worth noting that there is nothing particularly special about 
these modes to suggest they should be excited by physical processes (with the possible 
exception of Type IX, which may be written as maximal-wavelength gravitational waves 
on a closed background [19] )■ In particular, Type-V, VIIo and VII^ modes have 
symmetries which bear no relation to geodesies on the background manifold. While 
making the models interesting to CMB phenomenologists [TT], exactly this non- 
geodesic property makes them rather hard to motivate as candidates for "natural" initial 
conditions. On the other hand, the absence of any other way to transport anisotropic 
quantities around open spaces means that the Bianchi models must be taken seriously 
by any cosmologist who does not on physical grounds rule out the possibility of either 
an open or an anisotropic Universe. 

We closed by inspecting those Bianchi types which do not possess a strictly isotropic 
limit, concluding that all Bianchi types can appear isotropic and flat for an arbitrary 
period of physical time. We will presumably never be able to rule out observationally 
the possibility that we live in a Universe which is highly anisotropic on today's super- 
horizon scales. In particular, and in agreement with Wald's theorem [38j, the behaviour 
of all universes is seen to become increasingly isotropic if the late time is dominated by 
dark energy. 

The first-order evolution equations are extremely simple, and suitable for direct 
incorporation into CMB calculations; we have recently published temperature and 
polarization maps [12J by combining the results in the present work with our earlier 
derivation of the Boltzmann hierarchy for all nearly-isotropic Bianchi models [IT], 
excluding the types explained in Section [6] (the latter models would yield either a pure 
quadrupole or similar morphology to Type-V or Type- VIIo maps). Of course, although 
we tried to show (in Section |5| that our analysis does offer a different way of thinking 
about many previously known results, only a shadow of the complex behaviour which 
makes these models dynamically interesting [2] is captured in detail. 
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Appendix A. Alternative coordinate expressions for Killing fields and 
invariant vectors 

In Section |3j we gave expressions in terms of coordinate charts for the Killing fields of 
flat, open and closed maximally-symmetric spaces. We related these to the subset of 
fields which generate the Bianchi symmetries and listed the invariant fields for those 
cases. In this appendix we give, for completeness, these results in some alternative 
coordinate systems in the open and closed cases. As previously emphasized, these can 
be used to generate explicit expressions for the metric perturbations generated by each 
of the modes in Section HI 



Appendix A.l. Open space 



In order to express the results of Section 3.2 in terms of familiar coordinates, we define 
r, 0, by the following: 



coshr 

sinh r sin 9 cos < 



sinh r cos 9 
sinh r sin 9 sin i 



The metric induced by (12) is: 

ds 2 = dr 2 + sinh 2 r (dO 2 + sin 2 9 d(f) 2 ) 



(A.l) 



(A.2) 



With a little algebra, one may then obtain the KVFs in terms of the coordinate basis 



r ,e,. 



$,i = cos 9d r — coth r sin 9dg 

gVii h _ _y/}l cos 0Q r _|_ y/fa co th r sin Ode + d^ 

$,2 = sin cos 4>d r + (cos coth r — 1) cos dg + (cot — coth r esc 0) sin d<f, 

= sin sin 0<9 r + (cos coth r — 1) sin dg + (coth r esc — cot 0) cos 8$ (A. 3) 



The £2 and £3 in Type VII^ follow from equation ( |17| ). Near to r = 0, the Type-V KVFs 
reduce to Type I {i.e. translations along local Cartesian axes) as expected. Note also 
that the KVFs are well-defined everywhere despite the coordinate singularity at r = 0. 
The Type-V invariant vectors become 



,v 



cos cosh r — sinh r 



,v 



cosh 


r — 


cos sinh r 




sin 


COS 


cosh 


r — 


cos sinh r 




sin 


#sin0 


cosh 


r — 


cos sinh r 



8r 



8 r + 



8 r + 



csch r sin 



cosh r — cos sinh r 6 
cos (cos coth r — 1) 

cosh r — cos sinh r 
sin (cos # coth r — 1) 
cosh r — cos sinh r 



ft? — sin esc 6* csch rds 



dg + cos esc csch rc^ . (A. 4) 



The invariant basis for Type VII^ follows from equation (19) with x = — ln(coshr — 
sinh r cos 0). 
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An alternative form of the metric follows from foliating an open space into a series 
of flat 2-surfaces, resulting in 

ds 2 = dx 2 + e~ 2x (dy 2 + dz 2 ), (A.5) 

where x is as previously defined. The relation between this metric and the more usual 
cosmological metric is discussed in Ref. [35]. The coordinates are naturally adapted to 
the KVFs: starting at aP = 1, we reach any point on the open space by displacing a 
parameter distance x along £^ , followed by y along £^ and finally z along $^ . We then 
have 



.0 



1 



a" — -e x (y 2 + z 2 ) + coshx 

a 1 = ^e~ x (y 2 + z 2 ) + sinhx 
a 2 = e~ x y 

a 3 = e~ x z. (A.6) 

Because they are better adapted to the Bianchi symmetries, the KVFs and invariant 
fields look considerably simpler in these coordinates: 

e( = d x €Y = d x + yd y + zd z 

4 = e x d y £ = £ 2 VII 7^ = d y 



4 = e x d z $ = ff^/Vh = 3 Z (A.7) 



£vn h = y /h dx + ( yV £ _ z)dy + ( z ^jl + y ) dxi ( A . 8 ) 



where the ej llh can be obtained from (19) 



Appendix A. 2. Closed space 



Suitable coordinate expressions for the results in Section 3.3 may be expanded from 



equation (23). We include them here for the usual spherically-symmetric coordinate 
chart (for other common systems see |19j). 

a = cos r a 1 = sin r cos 9 

a 2 = sin r sin 9 cos a 3 = sin r sin 9 sin 0, (A. 9) 

for which the induced line element is 

ds 2 = dr 2 + sin 2 r (d9 2 + sin 2 9 d0 2 ) . (A. 10) 

The KVFs in this chart are 
2£j x = — cos 8d r + cot r sin 9d e + 

2£J, X = — cos sin 9d r — (cos 9 cos cot r + sin 0)<% + (cot r esc 9 sin — cos cot 9)8^ 
2£g X = — sin 9 sin 0<9 r + (cos — cos 9 cot r sin 0)<% — (cos cot r esc 9 + cot 9 sin 0)c^, 

(A.ll) 
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and the invariant fields are 

2e I 1 x = — cos 98 r + cot r sin 98 e — 8^ 

2e!, x = — cos sin 98 r + (sin — cos 9 cos 4> cot r)de + (cot r esc 9 sin + cos cot 
2e3 X = — sin 9 sin 0<9 r — (cos + cos 9 cot r sin 0)90 + (cot 9 sin — cos cot r esc 9)8^. 

(A.12) 



